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Abstract. According to recent observations, the existence of the dark energy has 
been considered. Even though we have obtained the constraint of the equation of 
the state for dark energy (p = wp) as —1 < w < —0.78 by combining WMAP data 
with other astronomical data, in order to pin down w, it is necessary to use other 
independent observational tools. For this purpose, we consider the w dependence of 
the non-Gaussianity of the density distribution generated by nonlinear dynamics. To 
extract the non-Gaussianity, we follow a semi-analytic approach based on Lagrangian 
linear perturbation theory, which provides an accurate value for the quasi-nonlinear 
region. From our results, the difference of the non-Gaussianity between w = —1 and 
w = —0.5 is about 4% while that between w = — 1 and w — —0.8 is about 0.9%. For 
the highly non-linear region, we estimate the difference by combining this perturbativc 
approach with N-body simulation executed for our previous paper. From this, we can 
expect the difference to be more enhanced in the low-z region, which suggests that 
the non-Gaussianity of the density distribution potentially plays an important role for 
extracting the information of dark energy. 
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1. Introduction 

According to recent observations for type la Supernovae pQ, the expansion of the 
Universe is accelerating. Combining measurements of Cosmic Microwave Background 
Radiation [2j and recent galaxy redshift survey jHj, researchers have concluded the 
Universe is almost flat, and we are forced to recognize the existence of a cosmo logical 
constant, or a kind of dark energy whose value is almost the same order of magnitude as 
the present density of the Universe [Ij. From the viewpoint of particle physics, however, 
it is quite difficult to explain such a tiny value that is 14 orders of magnitude smaller 
than the electroweak scale. The failure of theory to explain the present value of the 
cosmological constant is known as the "cosmological constant problem" j^J. 

In order to avoid this problem, many dark energy models have been proposed from 
various contexts. Roughly speaking, they can be classified into two approaches: one is 
to modify the gravitational law and the other is to introduce some exotic form of energy. 
For familiar examples, 1/R gravity models [0] and braneworld models [Zj belong to the 
former, while quintessence models 0IH], k-essence models [TU] and phantom models [TT] 
belong to the latter. 

Regardless of the attempts mentioned above, it is fair to say that a satisfactory 
explanation for dark energy has not yet been obtained. On the other hand, from the 
phenomenological viewpoint, it is important to constrain the effective equation of state 
of the dark energy w by observations. Even though we have obtained the constraint 
as w < —0.78 (95 % confidence limit assuming w > —1) 4J by combining WMAP 
data with other astronomical data, in order to pin down w, it is necessary to use other 
independent observational tools. 

For this purpose, we consider the statistical properties of the large-scale structure 
of the Universe based on the probability distribution function (PDF) of the cosmological 
density fluctuations. In the standard picture, the PDF of the primordial density 
fluctuations originated from quantum fluctuations that were stretched to large comoving 
scales during the inflation phase and are assumed to be a random Gaussian. It is, 
however, well known that even though PDF remains Gaussian as long as the density 
fluctuation is in the linear regime, it significantly deviates from the initial Gaussian 
shape once the non-linear stage is reached. 

Regarding the shape of the PDF, it has been shown that in standard cold dark 
matter (SCDM) models, the PDF is fairly approximated by lognormal distribution in 
a weakly non-linear regime while the lognormal PDF does not fit well in a highly 

non-linear regime [THj . After that, Plionis et al. ^3] and Borgani et al. analyzed 
the non-Gaussianity of the cluster distribution for several dark matter models including 
the low-density flat cold dark matter (ACDM) model by which the difference of non- 
Gaussianity between SCDM and ACDM models was clearly shown. (The comparison 
among them by N-body simulations for the quasi-nonlinear stage was also done in |16j.) 
In this paper, furthermore, we intend to investigate the dark energy model dependence 
for the non-Gaussianity of the PDF, in which we consider the simple constant-u> dark 
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energy models with w = —0.5, —0.8, —1.0, —1.2. 

In order to analyze the w-dependence of the non-Gaussianity of the PDF, we follow a 
semi-analytic approach based on Lagrangian linear perturbation theory for the evolution 
of the density fluctuation. This description provides relatively accurate values even in a 
quasi-linear regime in the structure formation scenario [121 H3 CHI UH1 1201 EH 1221 123 12H 
123 12H1 123 12H] • (Especially, for the case of the dust fluid, this approximation is called the 
Zel'dovich approximation (ZA)). Especially, in many suitable situations, it was shown 
that the Lagrangian approximation describes the evolution of density fluctuation better 
than the Eulerian approximation [23 E01 EH E2] • 

Another aim of this paper is the extension of our previous paper by including 
the dependence of the cosmic expansion of the background. This paper is organized as 
follows. In Sec. El we briefly present the evolution equations for the background fluid 
including the dark energy information. We then also present the perturbation equations 
based on Lagrangian linear perturbation theory (in Sec. El) and statistical quantities to 
investigate non-Gaussianity (in Sec. Hj). In Sec. El we provide the results and discuss the 
validity of our perturbative approach. Section El is devoted to conclusions. 

2. Background 

Here we briefly introduce the evolution equations for the background including the dark 
energy information. If we assume a homogeneous and isotropic Universe, the cosmic 
expansion law is described by the Friedmann equations and the energy conservation 
equation. 



where a is a scale factor of the Universe, H = a/ a is its Hubble parameter, JC is a 
curvature constant, and p and p are the total pressure and total energy of matter fields. 

As for matter fields, since we are interested in a stage much later than the radiation- 
matter equality, we consider matter fluid and dark energy, i.e., p = p m + Pde, where p m 
and pde are the matter fluid energy densities matter fluid and dark energy, respectively. 
In order to obtain the evolution of the energy densities for each component, we must 
specify the properties of dark energy which are strongly dependent on dark energy 
models. 

Even though there are some dark energy models in which dark energy couples 
explicitly to ordinal matter [Sj, dark energy interacts with ordinal matter only 
gravitationally in many models. We consider the case in which the energy conservation 
of each component holds independently. For these situations, since the pressure of the 
matter is negligible, its energy density evolves as usual: 




(1) 



(2) 
(3) 



p +3H(p + p) = 0, 



p m oca 



-3 



(4) 



Non-Gaussianity of the density distribution in accelerating universes 



4 



As for the evolution of dark energy, we can use the fact that it is possible to model 
dark energy as fluid in many dark energy models 5\. For this case, the dynamical 
properties of dark energy are determined through its effective equation of state w: 

Pde = w{pde)pde, (5) 

where, in general, w is an analytic function of the energy density, or equivalently, the 
scale factor. In principle, the evolution of energy density is determined by integrating 
the energy conservation equation: 

/ \ -3 ff dlna[w(a)+l] / r \ 

pDE[a) = pi,DEe Ja ' , (6) 

where p^de is an integration constant. The condition for the acceleration of the cosmic 
expansion is given as 

w < -i (7) 

where, if we are not to violate the strong energy condition, we should further impose 

-l<w. (8) 

However, since it is technically difficult to specify the time dependence of w by 
observations in near future, we concentrate on the case for constant w, in which case 
dark energy density evolves as 

p DE oc a- 3 ^ . (9) 

Assuming flat spacetime (JC = 0) and using Eqs. (J3J) and ©, we can rewrite the 
constraint Eq. (JTJ) as 

H 2 = Hi [Q m0 a- 3 + Q DEQ a' 3 ^} , (10) 

where we have defined the density parameters for the matter and the dark energy: 

o - 8nG n-n 

"(m, DE) = -^JpP(m, DE) ■ U-U 

The subscripts denote the corresponding quantities are evaluated at present. 
Under this constraint, we solve the Friedmann equation (Eq. (j2J|). 



3. Perturbation 

Next, we introduce the evolution equations for the density fluctuation of matter. Since 
we are interested in a scale much smaller than the horizon scale and non-relativistic 
matter fluid, we apply a Newtonian treatment. Notice that after this, p denotes the 
mass density, even though we used the same character as in the previous section. Here, 
we start with comoving coordinates along the cosmic expansion characterized by the 
solution a(t) of Eq. (fTOj) : 

r = a(t)x, (12) 

where r and x are physical coordinates and comoving coordinates, respectively. 
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In the comoving coordinates, the basic equations for cosmological fluid are described 
by the continuity equation, the Euler equation and the Poisson equation: 

^ +±V x -{v(l + 5)} = 0, (13) 
dv 1 , „ . a 1 _ 

— + - (v- V x )v + -v = -g, 14 
at a a a 

Vx9= -4irGp b a5, (15) 

where 5 = is the density fluctuation, p b is background density of the matter, v is 
the peculiar velocity, g = — ^V x $ is the peculiar acceleration and $ is the gravitational 
potential. Since we note non-relativistic fluid as matter, we can ignore the pressure and 
omit the suffix m. 

In this paper, for the perturbation, we adopt the Lagrangian picture rather than 
the Eulerian picture since we can extract the quasi-nonliniear nature of the structure 
formation even if we consider only the linear perturbation [2U EH EH] • For this purpose, 
it is necessary to define the comoving Lagrangian coordinates q in terms of the comoving 
Eulerian coordinates 

q = x + s(x,t), (16) 

where s is the displacement vector denoting the deviation from homogeneous 
distribution. While in Eulerian perturbation theory, the density fluctuation 5 is regarded 
as a perturbative quantity; in Lagrangian perturbation theory, the displacement vector 
s is regarded as a perturbative quantity, s can be decomposed to the longitudinal and 
the transverse modes: 

s = s L + s T , (17) 

Vxs L = 0, (18) 
V-s T =0, (19) 

where V means the Lagrangian spacial derivative. 

In Lagrangian coordinates, since we can solve the continuous Eq. ()13|) from Eq. (fTB^) 
exactly, we can obtain the density fluctuation as 

5 = J^ 1 - 1 , (20) 

where J is the determinant of the Jacobian of the mapping between q and r: dr/dq. 
It is worth noting that in Eq. (|2()jl . the density fluctuation is given in a formally exact 
form, even though we keep only the linear term in this paper. 

The peculiar velocity appearing in Eqs. (fT3*|) and (pPfj) can also be expressed in terms 
of the Lagrangian coordinate variable s as 

v = as. (21) 

The basic equations we shall solve are the linearized version of Eqs. ()14j) and 
(I15|) in which 5 and v are expressed with s from Eqs. (|2U|) and (}2Tj) . The solution of the 



Non-Gaussianity of the density distribution in accelerating universes 



6 



Lagrangian perturbation can be separated into a time-dependent part and a position- 
dependent part as: 

s^\t,q) = D(t)S^(q). (22) 

Since it is well known that there are no growing transverse mode solutions in the linear 
perturbation, in this paper we consider only the longitudinal mode in which the velocity 
field is irrotational. From the Kelvin circulation theorem, this is quite natural if it is 
generated only by the action of gravity 

The evolution equation for first-order solution is written as 

D + 2-D - AixGp b D = , (23) 
a 

while the spacial component is given by the initial condition. Zel'dovich derived a first- 
order solution of the longitudinal mode for dust fluid [T7j which provides the relation 
between the density fluctuation and the Lagrangian displacement as 

8(q) = -V-s L « . (24) 



It is worth noting that even though Eq. (j24j) is obtained by linear perturbation, the 
value of S remains accurate even in the quasi- nonlinear regime. 

For the low density flat Universe with dark energy whose equation of state is 
w = —1, i.e., 

VL m + Vt DE = 1, Vt DE = = const. , Pde = A , (25) 

the growing mode first-order solution can be expressed with the analytic function. 

h /5 2 
2 5l/h2 U'3 



D(t) = -B 1/h2 [-,-) , (26) 



h =421, (27) 
where B 1 / h 2 is an incomplete Beta function: 

B,(n,v)= f " p^il-pf-^p. (28) 



Jo 

For other dark energy models, we shall solve Eq. (|23J) numerically. 

In order to avoid the divergence of the density fluctuation in the limit of large k, 
however, just for a technical reason, it is necessary to consider the density field p(x; R) 
at the position x smoothed over the scale R, which is related to the unsmoothed density 
field p(x) as 



p(x;R) = I d i yW(\x-y\;R)p(y) 
d 3 k ~ 



W(kR)p(k)e- tk - x , (29) 
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where W denotes the window function and W and p represent the Fourier transforms 
of the corresponding quantities. In this paper, we adopt the top-hat window function, 

iy= 3(sinX - XCQSx) . (30) 

x 6 

Then, the density fluctuation 5(x; R) at the position x smoothed over the scale R 
can be constructed in terms of p(x; R). For simplicity we use 5 to denote 8(x; R) unless 
otherwise stated. 

4. Non-Gaussianity of the density fluctuation 

In order to analyze the statistics, we introduce a one-point probability distribution of 
the density fluctuation field P(S) (PDF of the density perturbation) which denotes the 
probability of obtaining the value 5. If 5 is a random Gaussian field, the PDF of the 
density perturbation is determined as 

P(5) = - e- s2/2 ° 2 

where a = ((8 — (5)) 2 ) and ( ) denotes the spacial average. 

If the PDF deviates from Gaussian distribution, in order to specify it, it is necessary 
to introduce the following higher-order statistical quantities [121 EH ESj : 

skewness : 7 - 1 1 
kurtosis : r\ 

which mean the display asymmetry and non-Gaussian degree of "peakiness," 
respectively. 

In Eulerian perturbation theory, the Gaussianity of the PDF completely conserves if 
we start with Gaussian distribution, because the density fluctuation is described by the 
product of the time and the spacial component and the spacial component never changes 
as long as we keep only linear terms. Once nonlinear terms are considered, however, 
its PDF deviates from the initial Gaussian shape because of the strong nonlinear mode 
coupling and the nonlocality of the gravitational dynamics. On the other hand, since 
in Lagrangian perturbation theory, the quasi-nonlinear information in the sense of the 
Eulerian picture can be extracted even by linear perturbation, we can expect to obtain 
nontrivial information for the skewness and the kurtosis after z ~ 5. See Fig. ^ 

5. Numerical results 




In this section, after mentioning the conditions we impose, we present our results. 

We consider here the Gaussian density field generated by COSMICS j^S], and we 
set the value of the density fluctuation and the peculiar velocity to be those obtained by 
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o I 1 1 1 1 1 

5 4 3 2 1 

Redshift z 

Figure 1. The dispersion of the density fluctuation in both Eulerian and Lagrangian 
linear perturbation for the model with w = —1. Here we choose the smoothing 
scale R = 8/i -1 Mpc for the top-hat window function. After z ~ 5, the difference 
of the dispersion appears between Eulerian and Lagrangian in which Eulerian linear 
perturbation no longer provides accurate values. Therefore we can expect to obtain 
the nontrivial information for the skewness and the kurtosis in this quasi-nonlinear 
region from Lagrangian perturbation as long as it provides accurate values. 

Lagrangian linear perturbation with the time evolution of the background to be a oc 
This is reasonable because at this time, the dark energy did not dominate the CDM. 
Accordingly, both the skewness and the kurtosis are less than 10~ 2 at initial (z ~ 30). 

For the computation of the Lagrangian perturbation, we set the parameters as 
follows: 

Number of grids : N = 128 3 , 

Box size : L = 128/i _1 Mpc (at a = 1) . 

Then we impose a periodic boundary condition. 

The cosmological parameters at the present time (a = 1) are given by j3] 



Q m =0.27, (31) 

Q DE = 0.73, (32) 

H = 71 [km/s/Mpc] , (33) 

a 8 = 0.84. (34) 



For the dark energy models, we set several equations of state. 

p = wp, (w = -0.5,-0.8,-1,-1.2), (35) 

where w = —1 corresponds to the cosmological constant. Even though w = —0.5 has 
already been excluded, we calculate this just for comparison. In the case of w — —1.2, 
as mentioned before, the strong energy condition is violated. 
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As for the smoothing scale R for the top-hat window function, we set 8/i _1 Mpc 
in the comoving Eulerian coordinates at the present time (a = 1). Because we apply 
Lagrangian linear perturbation, we must consider the validity of the perturbation. As 
we will show later, it is reasonable that we set the smoothing scale R = 8/i _1 Mpc if we 
are to discuss the structure formation with the perturbation. 

Based on 6 calculated from these conditions, we compute the statistical quantities, 
i.e., the dispersion, the skewness and the kurtosis at 51 time slices from z = 5 to z = 
whose time intervals are given as Az = 0.1. 

Before presenting the results, we must remind ourselves that since these calculations 
are based on Lagrangian linear perturbation theory, we can trust the results until 
caustics are produced, i.e., when a becomes ~ 1. 

Figure |2] shows the time evolution of the dispersion of the density fluctuation 
for several equations of state. The growth of the dispersion monotonously continues 
until z = 0. Because we choose a relatively large value for the smoothing length, the 
Lagrangian linear perturbation still has physical meaning. The tendency of structure 
formation can be discussed even if structure formation can no longer be fully described 
with linear approximation. 

We can identify the difference of the dispersion among several values of w. As the 
value of w becomes larger (approaches to 0), the dispersion becomes smaller. This can 
be explained as follows: For the model with larger w, dark energy dominates at earlier 
era, and the expansion of the Universe starts to accelerate at earlier era, while it can 
be shown that the growth of the density fluctuation is smaller in the accelerating stage 
than in the matter-dominant stage. 

Here we consider high-z region, i.e., z > 2. In our previous paper jSH], we 
compared the evolution of the dispersion between N-body simulation and Lagrangian 
linear perturbation. There, the difference of the dispersion, the skewness, and the 
kurtosis between N-body simulation and the Lagrangian linear perturbation stayed small 
until o ~ 0.3. From the past analyses, the Lagrangian linear perturbation describes 
structure formation rather well until z > 2, where the dispersion is a ~ 0.3. 

Withrespect to non-Gaussianity, Figs. HU and |U show the evolution of the skewness 
and the kurtosis, respectively. Like the dispersion, the growth of the skewness and the 
kurtosis continues monotonously. 

For the reason mentioned above, our computations seem to be reliable until z ~ 2. 
For the case where w = —0.5, the value of both the skewness and the kurtosis is obviously 
less than that for the cases where w = —1.2, w = —1.0 and w = —0.8 (about 4% and 
about 8%, respectively). On the other hand, the differences of both the skewness and 
the kurtosis among w = —1.2, w = —1 and w = —0.8 are less than 2% at z = 2. Even 
though these are very small values, they are of almost the same order as the relative 
ratio of the dispersions. Therefore it can be expected that non-Gaussianity can play a 
potentially important role in distinguishing the value of w from the observation of high-z 
galaxies. At that stage, the fitting functions for j w (z), r] w (z) for the quasi-nonlinear stage 
are useful. Since it seems that the growth of the dispersion and the non-Gaussianity are 
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Figure 2. The dispersion of the density fluctuation in several dark energy models. 
The growth of the dispersion continues monotonously. The difference of the dispersion 
between w = — 1 and w — —0.8 is about 0.9% at z — 2 until which we regard Lagrangian 
linear perturbation as valid. In the same way, the difference of the dispersion between 
w = — 1 and w = —1.2 is about 0.4% at z = 2. 

2.5 



2 
1.5 
1 

0.5 



5 4 3 2 1 

Redshift z 

Figure 3. The skewness of the density fluctuation in several dark energy models. The 
growth of the skewness continues monotonously. Although we can easily distinguish 
the case where w — —0.5 from the results, the difference of the skewness between 
w = —1 and w — —0.8 is about 0.9% at z = 2 until which we regard Lagrangian linear 
perturbation as valid. In the same way, the difference of the skewness between w — — 1 
and w = —1.2 is about 0.4% at z = 2. 

related, it is worth trying to find the functions based on our results. 

Finally, if we want to know precise information for non-Gaussianity after z ~ 2, 
we must consider nonlinear dynamics by such as N-body simulations. Even though we 
do not compute them in this paper, judging from Fig.l in our previous paper, jHH], in 
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w=-0.8 
w=-0.5 
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Figure 4. The kurtosis of the density fluctuation in several dark energy models. 
The growth of the kurtosis monotonously continues. As well as the distinction of the 
skewness, the difference of the kurtosis between w = — 1 and w = —0.8 is about 1.8% 
at z = 2 until which we regard the Lagrangian linear perturbation is valid. In the 
same way, the difference of the kurtosis between w — — 1 and w = —1.2 is about 0.9% 
at z = 2. 

which we compare the results of Lagrangian perturbation and N-body simulation in the 
Einstein-de Sitter Universe, we can expect the actual non-Gaussianity is strengthened by 
a order of magnitude. If this happens also for accelerating Universes, it becomes easier 
to distinguish w from the observations of the near galaxies, which is worth investigating. 

6. Summary 

To clarify the nature of the dark energy is, without doubt, one of the most important 
tasks in modern cosmology. From the phenomenological viewpoint, it is important 
to constrain the effective equation of state of dark energy w by observations. Even 
though we have obtained the constraint as w < —0.78 (95 % confidence limit assuming 
w > —1) [1] by combining WMAP data with other astronomical data, to pin down w, 
it is necessary to use other independent observational tools. 

For this purpose, we consider the non-Gaussianity of the density fluctuation which 
is generated by the non-linear dynamics. In this paper, we follow a semi-analytic 
approach based on Lagrangian linear perturbation theory for the evolution of the density 
fluctuation. In this theory, we can extract the quasi-nonlinear feature of the density 
fluctuation even though we keep only linear terms. In terms of the density fluctuation 
we obtain, we compute the skewness and the kurtosis which are statistical quantities 
denoting the degree of non-Gaussianity as well as the dispersion, which is the only 
parameter for Gaussian distribution. By considering several constant- 1 ?/; dark energy 
models, we present the w dependence of non-Gaussianity. Because of the validity of the 
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perturbative approach, we regard the results until z ~ 2 as reliable. 

According to our calculation, at fixed time (z), the dispersion becomes smaller as w 
becomes larger (approaches to 0). This is because the larger w is, the earlier the time at 
which the expansion of the Universe starts to accelerate and the growth of the density 
fluctuation is smaller than in the matter-dominant Universe. 

The skewness and the kurtosis also become smaller as w becomes larger (approaches 
to 0), which suggests they have some correlation with the dispersion. The difference of 
non-Gaussianity between w = —1 and w = —0.5 is obvious (about 4%) for 2 < z < 5, 
while it is hard for w to distinguish —1.2, —1 or —0.8. The differences of both the 
skewness and the kurtosis between w — — 1 and w = —0.8 as well as w = —1.2 and 
w = —1.0 are quite small (about 0.9%), while the correspondent differences of the 
dispersion are the almost same as these. To specify the value of w, the fitting functions 
for j w (z), T] w {z) for the quasi-nonlinear stage which can be constructed from our results 
are useful. 

For regions nearer than z — 2, it is necessary to pick up non-linear information 
for obtaining the degree of non-Gaussianity. Even though we have not computed these 
in this paper, from our previous results comparing N-body simulation with Lagrangian 
perturbation theory, it can be expected that the actual value of non-Gaussianity becomes 
larger than the value obtained by Lagrangian linear perturbation. For this case, it 
becomes easier to distinguish the value of w from the observation of near galaxies, 
which is worth continuing to consider. 

It is necessary to mention the possibility of examining our results. Recently, several 
galaxy redshift survey projects have been progressing [371 EE1 EH] . In these projects, 
many galaxies within a region (z < 0.3 for 2dF, z < 0.5 for SDSS, and 0.7 < z < 1.4 for 
DEEP2) have been observed. These projects show the latest distribution of galaxies, 
which form strongly nonlinear structure. For the next generation of spectroscopic 
surveys, many high-z galaxies will be observed. For example, the Cosmic Inflation 
Probe (CIP) project [10] is planned to detect objects between 3 < z < 6.5. Another 
project, the Kilo- Aperture Optical Spectrograph (KAOS), has also been proposed 
One of the primary scientific purposes of KAOS is the determination of the equation of 
state of dark energy by objects until z < 3. Another purpose of KAOS is to observe 
the growth of structure. The projects leaders propose that the linear growth factor of 
matter density fluctuations can also provide sensitivity to dark energy properties. 

The skewness and the kurtosis of the distribution of galaxies in the SDSS Early 
Data Release had been computed [12] ■ The results suggested the SDSS imaging data 
can enable us to determine the skewness and the kurtosis up to 1% and less than 10%. 
Therefore in future projects, we can expect to determine the skewness and the kurtosis 
with high accuracy. From our results, the kurtosis of the density fluctuation seems more 
sensitive for w than the dispersion of the fluctuation. In other words, although the 
dispersion reflects the linear growth factor, the kurtosis shows the nature of nonlinear 
growth in addition to the linear growth. Therefore, structual non-Gaussianity of would 
show more detailed information on the nature of dark energy. It is worth noting that 



Non-Gaussianity of the density distribution in accelerating universes 



13 



for such a high- z region, the results obtained in this paper based on our perturbative 
approach will play an important role. 

Finally, it is worth noting the primordial non-Gaussianity which is generated by 
inflation. From WMAP observations, a limit for the non-linear coupling parameter 
has been established jl3] which does not deny the existence of non-Gaussianity in the 
primordial density fluctuation. Even though we cannot disentangle them generally, we 
can calculate and compare non-Gaussianity of large-scale structure from both Gaussian 
and non-Gaussian initial conditions. This requires further investigation, and we hope 
to report results in a separate publication. 
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